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A kinetic Monte Carlo method was used to simulate the diffusion of reptating polymer chains 
across the interface. A time-resolved fluorescence technique conjunction with direct energy trans- 
fer method was used to measure the extend of diffusion of dye labeled reptating polymer chains. 
The diffusion of donor and acceptor labeled polymer chains between adjacent compartments was 
randomly generated. The fluorescence decay profiles of donor molecules were simulated at several 
diffusion steps to produce mixing of the polymer chains. Mixing ratios of donor and acceptor la- 
beled polymer chains in compartments were measured at various stages (snapshots) of diffusion. It 
was observed that for a given molecular weight, the average interpenetration contour length was 
found to be proportional to the mixing ratio. Monte Carlo analysis showed that curvilinear diffusion 
coefficient is inversely proportional to the weight of polymer chains during diffusion. 



I. INTRODUCTION 



The diffusion of polymer chains across polymer- 
polymer interfaces has been of interest for more than a 
decade. li^J One of the reason for the interest is that poly- 
mer diffusion across an interface is found to be important 
in technological processes, such as sintering of polymer 
powders, development of the strength of polymer powder 
by compression molding and annealing and the formation 
of latex films. Latex- film formation has been considered 
in the literature for over 50 years and is frequently used 
in the modern industry.^ JJ The process of latex film 
formation has been divided into several stages. The gen- 
erally accepted mechanisms consist of: (i) evaporation 
which brings the particles into some form of close pack- 
ing; (ii) deformation of particles which leads to a struc- 
ture without voids, although with the original particles 
still distinguishable; and (iii) diffusion of polymer chains 
across particle-particle boundaries, yielding a continuous 
film with mechanical integrity. Voyutskii proposed that 
physical contact between latex particles would not pro- 
duce a mechanically strong continuous film if no external 
effect is applied.^ In other words, in order to obtain a 
stable film it is necessary that the segments of polymer 
chains diffuse from one particle to another by forming 
a strong linkage between them. He described the pro- 
cess of diffusion with the word of "autohesion" . For sev- 
eral decades, after Voyutskii's paper there have been only 
speculations about the process of diffusion during latex 
film formation. The developments in neutron scattering 
techniques have enabled the carrying out of experiments 
to study diffusion of polymer molecules across particle- 
particle boundaries. Important progress was made when 
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small angle neutron scattering (SANS) technique was ap- 
plied to study diffusion between particles of hydrogenated 
and deuterated acrylic latex particles. IS Later, it was re- 
ported that increasing molecular weight and incompati- 
bility lower the diffusion rate J 7 ' SANS was employed to 
measure the extend of diffusion in polystyrene (PS) la- 
tex having high molecular weight and small particle size, 
where an increase in the radius of gyration (R g ) of the 
polymer was observed when the system was heated above 

itS TgM 

The non-radiative direct energy transfer (DET) 
method conjunction with time resolved fluorescence 
(TRF) tecnique was first used to study polymer diffusion 
across particle-particle boundaries^ to monitor concen- 
tration profiles of donor and acceptor dyes attached to 
polymers that are located initially in separate particles. 
As polymer diffusion occurs, mixing of donor (D) and 
acceptor (A) dyes can be measured by an increase in en- 
ergy transfer between them. Early measurements using 
DET on poly (methyl methacrylate) (PMMA) latex par- 
ticles prepared by nonaqueous dispersion polymerization 
found diffusion coefficients of the order of 10~ 15 cm 2 s~ 1 
at temperatures between 400 — 450 K. The film formation 
of a poly (butyl methacrylate) latex prepared via emul- 
sion polymerization in water, was studied using the same 
technique. The diffusion coefficients (10 _16 cm 2 s _1 ) were 
determined using spherical diffusion model, and found to 
be dependent on both time and temperature. A decrease 
in diffusion rate with time was attributed to the effect 
of low molecular weight chains near the particle-particle 
surface dominating at early times. The DET method 
of data analysis was later developed for melt pressed 
PMMA particles. li2J It is observed that mass transfer in- 
creased with time to a power of 0.5 as in Fickian diffu- 
sion model for low molecular weight polymers. For high 
molecular weights, there is a distinct 0.25 power depen- 
dence that cannot be explained by a Fickian or reptation 
model. 

Two reviewsIif^AJ summarized experimental ap- 



proaches and results in studying diffusion in latex sys- 
tems up to the year of f 993 and 1994, respectively. Both 
papers reviewed the evidence of diffusion of polymers 
at particle-particle junction using transmission electron 
microscopy, SANS, atomic force microscopy and fluores- 
cence techniques. Data analysis using DET and a flu- 
orescence technique was later improved by taking into 
account the donor and acceptor concentration profiles 
during polymer diffusion, where a uniform acceptor con- 
centration was considered around a donor in thin slices 
or shells. Ii2f22) Further a model for DET was developed 
which considered the heterogeneity in the donor and ac- 
ceptor concentration profiles,!^ where the diffusion coef- 
ficient of polymer chains obtained by different DET mod- 
els. DET studies in latex blends, where one phase is far 
below its T g and does not undergo any significant dif- 
fusion, found that the magnitude of energy transfer is 
proportional to the interfacial area. It was observed that 
organic solvents, which plasticizc the latex, enhance the 
rates of polymer diffusion and the diffusion rate increased 
in the presence of cosurfactants.tiS At the times longer 
than the tube renewal time (escape time), the activation 
energy for diffusion was found 30 kcal/mol in the pres- 
ence of cosurfactants, which is 48 kcal/mol in the neat 
polystyrene. Monte Carlo studies on diffusion of donor 
and acceptor dyes between the adjacent compartments 
were modeled by using Brownian motion. li&2£22J 

In this work Kinetic Monte Carlo (KMC) method was 
used to simulate the diffusion of polymer chains between 
the adjacent compartments of a cube where free energies 
and the potentials between molecules were not consid- 
ered. These chains are labeled with either donor or ac- 
ceptor dyes. Reptation of polymer chains were simulated 
using KMC method which produced more realistic diffu- 
sion results compared to previous studies. life 22) The 
diffusion of donor and acceptor labeled polymer chains 
between adjacent compartments was randomly gener- 
ated where each chain reptates according to de Gennes 
model. The decay of the donor fluorescence intensity, 
I(t) by DET was simulated at several diffusion steps and 
a gaussian noise was added to generate the time resolved 
fluorescence decay data. I(t) decays were then fitted to 
the phcnomcnological equation 



TU) - B lC M-t/r - C(-) 1 / 2 ) + B 2 eM-t/r ) (1) 



1(0) 



TO 



where To is the fluorescence lifetime of donor and B\ and 
i?2 are the pre-exponential factors. Here it is useful to 
define the mixing ratio K representing the order of mix- 
ing during diffusion of the donor and the acceptor labeled 
polymer chains as 
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K was measured in terms of B\ and £? 2 for various molec- 
ular weights, M of polymer chains. Results were in- 



FIG. 1: A cartoon representation of a reptating polymer chain 
in a tube at different times (a,b,c and d) . The memory of the 
initial tube is lost in Figure Id. The picture also demonstrates 
the motion of minor chain during reptation. 



terpreted in terms of average interpenetration contour 
length, £(t r ) and curvilinear diffusion coefficient, D\. 



II. THEORETICAL CONSIDERATIONS 

A. Molecular Motion of Polymer Chains 

The motion of individual chains in bulk polymeric ma- 
terials or concentrated solutions of linear random-coil 
(Gaussian) chains has been modeled by the reptation 
theory of de GennesIiS) and Edwards. 122) In this model 
the polymer chain is confined to a tube which presents 
topological constraints to lateral motion of monomers im- 
posed by the neighboring chains via entanglements. The 
motion of the chain is restricted to the curvilinear length 
of the tube, which is shown in Figure 1. As shown in 
Figure la, at time zero, the chain is in its initial tube 
which performs Brownian back and forth motion. Since 
the chain ends are free to move in any direction away 
from the tube, the memory of the initial tube position 
in space is gradually lost as shown in Figures lb and lc. 
Figure Id presents the final stage of the reptation, i.e. 
just before the total escape of the chain from its original 
tube. At the time T r , the escape time, the chain escapes 
or forgets its original configuration. As the chain is in a 
bulk polymeric system, it always creates a new tube af- 
ter escaping from the previous one. The length of chain, 
£(t r ) which escapes from the initial tube is also a random- 
coil chain which obeys Gaussian statistics and called the 
minor chain. I2if22) During the process of reptation the 
length of the minor chain grows which is represented in 
Figure 1 by the growing spherical envelope. The mean 
square escape length of the minor chains, < £ 2 > is cal- 
culated by the following relationl2i) 
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FIG. 2: Schematic picture of the bond and dihedral angles for 
the polymer chains. The bond angle <j> is equal to 60°, and 
the dihedral angles, namely 9, can be randomly chosen. 



< £ 2 >= Wtt^DiU . (3) 

where t r is the diffusion time and D\ is the curvilin- 
ear one-dimensional diffusion coefficient. The molecular 
weight dependence of Di and T r are given by the follow- 
ing relations 

Di - M~ l (4) 

and 

T r ~ M 3 (5) 

respectively. 

The average interpenetration contour length, £(t r ) of 
chain segments which have diffused across the interface 
is obtained from the minor chain model aslSSJ 

£(t r ) ~ M- 1 ' 2 ^ 2 (6) 



too ~ M (7) 

This property has the same scaling relation as the minor 
chain length, < I >. 

B. Direct Energy Transfer Method 

Time resolved fluorescence (TRF) in conjunction with 
DET method monitors the extent of diffusion of donor 
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FIG. 3: The plot of a) escape time, T T against the number 
of reptation steps, d, b) (Tr) 1 / 3 versus molecular weight, M, 
for different reptation steps, d, for a polymer chain containing 
200 monomer units. 



(D) and the acceptor (A) labeled polymer molecules. If 
the sample is made of a mixture of D and A labeled poly- 
mer chains where the diffusion takes place, after a period 
of time when the donor fluorescence profiles are mea- 
sured, each decay trace provides a snapshot of the extent 
of diffusion. 12) This sample is considered to be composed 
of three regions; unmixed D, unmixed A and the mixed D 
- A region. The above model was first empirically intro- 
duced by the two component donor fluorescence decay. 123) 
When donor dyes are excited using a very narrow pulse 
of light, the excited donor returns to the ground state ei- 
ther by emitting a fluorescence photon or through the 
nonradiative mechanism. For a well behaved system, af- 
ter exposing the donors to a short pulse of light, the fluo- 
rescence intensity decays exponentially with time. How- 
ever, if acceptors are present in the vicinity of the excited 
donor, then there is a possibility of DET from the ex- 
cited donor to the ground state acceptors. In the classical 
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FIG. 4: The plot of a) mean square escape length, < £ 2 > of 
minor chain versus diffusion time, t r , b) The inverse slope of 
< I 2 > versus molecular weight, M. For a polymer chain of 
200 monomer units. 



problem of DET, neglecting back transfer, the probabil- 
ity of the decay of the donor at r k due to the presence of 
an acceptor at n is given byl22J 
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Here nu and ua represent the distribution functions of 
donors and acceptors. In the thermodynamic limit equa- 
tion (|10|) becomes 
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— = exp(-t/T ) — J n D (r k )dr k 

x exp(- J n A (n)dri(l - exp (-w ik t))) .(11) 

This equation can be used to generate donor decay pro- 
files by Monte-Carlo techniques. It is shown that the 
equation i|ll|) reduces to a more simple form which can 
be compared to the experimental data.t^J The simplifi- 
cation is summarized below for clarity. Changing to the 
coordinate = j-j — r k leads to, 
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where R g is an arbitrary upper limit. Placing a particular 
donor at the origin and assuming that the mixed and 
unmixed regions are created during diffusion of D and A, 
equation Ijl2|) becomes 
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where 



P k (t) = exp[-t/r - Wikt] 



(8) 



where Wi k is the rate of energy transfer given by Forster 
as 
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Here i?o represents the critical Forster distance and k 
is a dimensionless parameter related to the geometry of 
interacting dipole. If the system contains Np donors and 
N A acceptors, then the donor fluorescence intensity decay 
can be derived from the equation @ and given 

by [io, nj 
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represent the fraction of donors in mixed and unmixed 
regions respectively. The integral in equation H13|l pro- 
duces a Forster type of function!^ 211 



Na 1 f R » t 
TT-— / n A (r ik )dr lk exp(-w ik t) = exp(-C( — ) 1/2 ) 

t=i Na Jo t ° 

(15) 

where C is proportional to acceptor concentration. Even- 
tually, one gets equation Q for the fluorescence intensity. 
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FIG. 5: Snapshots of diffusion between adjacent compart- 
ments of the cube with the size £xix L(L — 30). Donor, D 
labeled chains in the right compartment are presented in black 
for clarity. The chains consist of 10 units and the number of 
chains in each compartment is taken to be 20. The grey col- 
ored chains represent the ones that are labeled by acceptors, 
A. 



FIG. 6: a) Donor decay profiles at several diffusion steps for 
M = 100 and d — 10. Where the mixing ratios are taken as 
I) K = 0.1, II) K = 0.5, III) K = 0.9, b) Noisy decay profiles 
for the above picture. 



III. RESULTS AND DISCUSSION 

A. Freely Rotated Chain Model for Reptation 

The polymer chains were constructed from the rod like 
segments which were connected to each other with a bond 
angle, <f> = 60°. The dihedral angle, 9 was chosen ran- 
domly as seen in Figure 2. The length of a segment was 
taken to be unity and the molecular weight, M of the 
chain was assumed to be proportional to the number of 
rods (monomers) that chain contains. Monte Carlo simu- 
lations were performed so that the chains move according 
to the reptation model, where for each reptation step, 
while a segment from the tail disappears, another seg- 
ment is added to head or vice versa. In other words the 
chain is being trapped in a hypothetical tube, leaves the 
tube in a randomly chosen directions, creating a segment 
of minor chain as shown in Figure 1. A unit of diffu- 
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FIG. 7: a) Variation of K with respect to diffusion time, t r 
for different molecular weights, I) M = 100, II) M = 150, III) 
M = 200, IV) M = 250 for d = 15, b) The plots K versus 
diffusion time, t r for different d values; I) d = 25, II) d — 20, 
III) d = 15, IV) d = 10, V) d = 5 and VI) d=l, for M = 200. 



FIG. 8: a) K versus t r 1//2 plots for different values of M, b) 
K versus t r 1//2 plots for different values of d. Slopes of linear 
fits produce a in equation 1171 . 



sion time, t r was defined for a single reptation step of a 
chain that has 100 monomers. For instance single repta- 
tion step of a chain with 200 monomers takes 2 units of 
t r . The number of reptation steps in chosen direction for 
each simulation was defined as d which controls the es- 
cape time i.e. an increase in d shortens the escape time, 
T r . Figure 3a presents the behavior of escape time, T r 
against d for a polymer chain which has 200 monomer 
units. Here T r was determined with a separate program 
which uses the same reptation algorithm as in the sim- 
ulation programs. In this program a single chain rep- 
tates until it escapes from its initial tube and then T r is 
recorded. The average T r is calculated from several runs. 
It is seen in Figure 3a that as d is increased T r decreases 
i.e. chain escapes quicker from its initial tube. The re- 
lation between T r and molecular weight, M for different 
d values are shown in Figure 3b, where linear behavior 

1/3 

of T r versus M suggest that polymer chains reptate in 
accordance with equation (J5J in our KMC simulations. 



The steepest curve for d = 1 predicts that corresponding 
polymer chain needs longer time than others to escape 
from its tube. However when chain reptates with d — 25 
escape time, T r is very short for all molecular weights. 

The minor chain growth was monitored by observ- 
ing the mean square escape length < £ 2 > versus dif- 
fusion time, t r where < £ 2 > was determined using the 
same separate program which determines T r . This pro- 
gram reptates a single chain and records the minor chain 
length. The average was calculated after several runs of 
this program. Results are shown in Figure 4a for a chain 
contains 200 monomer units, where it is seen that the 
region between t r — to t r = 0.3T r presents linear rela- 
tion which accords with equation ©. The slope of the 
curves in the linear region is proportional to the curvilin- 
ear one dimensional diffusion coefficient, D\ which were 
produced for various molecular weights, M. The inverse 
of the slopes, (slope) -1 are plotted versus M in Figure 4b 
where the perfect linear relation is in accord with equa- 
tion gj) i.e. Di - M _1 . 
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B. Energy Transfer for Chain Diffusion 

Diffusion of D and A labeled chains between adja- 
cent compartments was simulated using above KMC al- 
gorithm. The sample cube with the side, L = 500 units, 
was divided into two equal compartments which contain 
donor, D and acceptor, A labeled chains respectively, 
presented in Figure 5. Each compartment contains 500 
chains and one percent of monomers in each chain labeled 
with donors or acceptors. Here one has to be noticed that 
chains are chosen quite short and the system is very di- 
lute for the reptation model. In fact, here the reptation 
algorithm is designed so that it operates as if the chains 
are in the melt system. This algorithm is chosen to avoid 
the immediate decay of donors where the chains are kept 
far away from each other and still perform the reptation 
motion. One may attemp to solve this problem by filling 
up the system with unlabeled chains by designing more 
realistic algorithm, which, however then spends tremen- 
dous computer time to perform reptation motion for all 
chains. Since these unlabeled chains have no contribution 
to the DET, using the reptation algorithm in the dilute 
system saves considerable amount of computer time and 
still mimics the realistic reptation motion. The measured 
reptation parameters have shown that the chosen dilute 
system and the reptation algorithm work quite well and 
the attempt to study the DET during diffusion is reliable. 

The decay of donor intensity by DET was simulated for 
the possible configurations at the end of each 10 4 steps 
of diffusion, therefore the snapshots of diffusion processes 
can be monitored quite clearly and accurately. Reflected 
boundary conditions were used from the sides of the cube. 
For the configurations of d = 1 the decay of donor inten- 
sity was simulated after 0.5 x 10 5 — 10 5 steps of diffusion 
because diffusion process for this case takes much longer 
time. Snapshots of the diffusing polymer chains between 
adjacent compartments in a cube are shown in Figures 
5a, b and c at various diffusion steps. The D and A la- 
beled chains are colored in black and grey, respectively, 
for clarity. 

The donor decay profiles were generated using equa- 
tion (|lip. The Wik values for each donor-acceptor pair 
were obtained from equation © using a Forster distance 
of 26 units. The parameter k 2 was chosen as 0.476, a 
value appropriate for immobile dyes,ti3 and the donor 
lifetime To was taken as 44ns. Then equation was 
used to simulate the donor intensity I(t) i.e. fluores- 
cence decay profile. 1(0) = 2 x 10 4 was chosen and the 
decay profiles were obtained for a 250ns interval, divided 
into 250 channels of Ins each. The fluorescence decay 
profiles I(t), for the polymer chains with M = 100 and 
d = 10 are given in Figure 6a at various diffusion steps, 
i.e. various mixing ratios, K. In order to obtain more 
realistic decay profiles, gaussian noise was added to the 
original decay profiles using Box, Muller and Marsaglia 
(1958) algorithm J 28 ' Here, one may also take into ac- 
count the effect of the lamp profile when calculating the 
decay profiles. To do so the decay profiles gener- 
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FIG. 9: a) The dependence of a on molecular weights, M, for 
various values of d, b) Plots of the variation of a with d, for 
different values of M. The linear relation between a and d 
are shown by the regression lines. 



ated by the Monte Carlo simulation should be convolved 
with experimental lamp profile, then the experimentally 
measured ip(t) is obtained by convolution of I(t) with the 
instrument response function L(t), as 



L(t)I(t - s)ds 



(16) 



In this work generated decay 
I(t)J 16 ' This assumption is 
r5 function light source (e.g. 
lamp profile. In this case no 
equation H16|) produces I(t). 
of donors for a S function lig 
Figure 6b at various diffusion 



profiles were used, namely 
valid if one uses a delta, 

a very fast laser) as the 
convolution is needed and 

The noisy decay profiles 
lit source are presented in 
steps, similar to Figure 6a. 



Extend of Diffusion 



In order to calculate the mixing ratio, K defined in 
equation @ , the decay profiles were fitted to equation 
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FIG. 10: a versus M" 1/2 plots for different values of d. The 
linear regression lines show the goodness of fit to equation 
Gg- 



using Levenberg-Marquartt22J algorithm. During curve 
fitting process C(=l) and To (=44) were kept constant 
and B\ and Bi values were varied. Since C and To were 
fixed, the fitting procedure directly produced B\ and B2 
values, which are a pre-exponential factor in mixed and 
unmixed regions. At early times of diffusion B2 domi- 
nates by presenting low mixing, however at later times 
B\ increases and dominates the mixing ratio. Here ba- 
sically simulation of decay curves are essential in calcu- 
lation of B\ and Bi values. More than 10 4 decay pro- 
files were fitted and the goodness of fitting was accepted 
as x 2 < 1-5- The B\ and B 2 values were used to ob- 
tain K values at diffusion steps during reptation of poly- 
mer chains. Figure 7a presents the variation of K with 
respect to the diffusion time, t r for chains at different 
molecular weights, M (100, 150, 200, 250) for d = 15. As 
seen in Figure 7a mixing of diffusing chains are much 
faster for the low molecular weight samples (I, II) than 
high molecular weight samples (III, IV). When the molec- 
ular weight is increased diffusion slows down as expected. 
Diffusion for different d values is also shown in Figure 7b, 
K versus diffusion time are plotted for various d values 
(1,5,10,15,20,25) for chains at M = 200, where it is 
seen that increase in d cause increase in mixing of dif- 
fusing chains. In order to quantify the above results the 
following equation was employed 



Here we intent to elaborate the equation ©, where the 
molecular weight, M can be related with a in equation 
(|T7|) . The fits of the data in Figures 7a and 7b to equation 
(|17|l are given in Figures 8a and 8b, respectively. The 
slopes of the linear relations in Figure 8 produce "a" 
values which are plotted versus M and d in Figures 9a 
and 9b, respectively. Curves in Figure 9a predict that as 
M is increased, a decreases. In order to determine the 
behavior of a with respect to M; an a versus M -1 / 2 plot 
is presented in Figure 10 where the following relation is 
obeyed. 



M~ 



-1/2 



(18) 



Here one may speculate that mixing ratio, K is most 
probably proportional to the average interpenetration 
contour length, l(t r ) according to equation @. In other 
words K is the measure of £(t r ) which is quite important 
in determining the mechanical properties of polymeric 
materials. 122 

On the other hand the plot in Figure 9b suggests that 
as d is increased, a increases which predicts that as the 
chains reptate faster, the mixing ratio, K increases as 
expected. It is evident that since small chains reptate 
faster, they can mixe much quicker than high molecular 
weight chains. Here one may also argue that the root 
square of equation suggests that D\ ~ A/ -1 relation 
holds for our simulation results. 



K = at 



1/2 



(17) 



IV. CONCLUSION 

In conclusion this paper has presented a Kinetic Monte 
Carlo method which can be used to simulate the inter- 
diffusion of reptating polymer chains across the inter- 
face. It was shown that Equation (1) can be succesfully 
applied to measure the extend of interdiffusion of reptat- 
ing donor and acceptor labeled polymer chains. Monte 
Carlo results have shown that the curvilinear diffusion 
coefficient is inversely proportional to the weight of the 
polymer chains and for a given molecular weight the av- 
erage interpenetration contour length was found to be 
proportional to the extent of interdiffusion. 
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